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ABSTRACT 

We explore the idea that the observed variations in the peak luminosities of Type la supernovae originate in 
part from a scatter in metallicity of the main-sequence stars that become white dwarfs. Previous, numerical, 
studies have not self-consistently explored metallicities greater than solar. One-dimensional Chandrasekhar 
mass models of SNe la produce most of their 56 Ni in a burn to nuclear statistical equilibrium between the mass 
shells 0.2 Mq and 0.8M Q , for which the electron to nucleon ratio Y e is constant during the burn. We show 
analytically that, under these conditions, charge and mass conservation constrain the mass of 56 Ni produced to 
depend linearly on the original metallicity of the white dwarf progenitor. Detailed post-processing of W7-like 
models confirms this linear dependence. The effect that we have identified is most evident at metallicities larger 
than solar, and is in agreement with previous self-consistent calculations over the metallicity range common 
to both calculations. The observed scatter in the metallicity (^Zq-3Z ) of the solar neighborhood is enough 
to induce a 25% variation in the mass of 56 Ni ejected by Type la supernova. This is sufficient to vary the 
peak V-band brightness by |AMy| « 0.2. This scatter in metallicity is present out to the limiting redshifts of 
current observations (z < 1). Sedimentation of 22 Ne can possibly amplify the variation in 56 Ni mass to < 50%. 
Further numerical studies can determine if other metallicity-induced effects, such as a change in the mass of 
the 56 Ni-producing region, offset or enhance the variation we identify. 

Subject headings: nuclear reactions, nucleosynthesis, abundances — stars: abundances — supernovae: general 
— galaxies: abundances — galaxies: evolution 



1. INTRODUCTION 

The maximum luminosity of Type la supernovae (SNe la) is 
a key ingredient of their analysis and is also essential for their 
use as distance indicators in c osmology {Fillipenko, Il997t 
iBranchl 19981 lLeibundgu tl2001l) . For the nearby SNe la with 
Cepheid-determined distances, the overall dispersion in the 
peak magni tude measurements is rather small, about 0.5 mag- 
in B and V JFillipenkoll997tlSaha etafl i ll 9991) : lGibson et alJ 
2000). When the sample is enlarged to include more distant 
SNe la, there are several sublumin ous events that broaden the 
variation to about 1 mag, in B jH^mjLy_£Lailll996; Riess et al. 
ll998t lPerlmutte r et alJl999HPhillips et alJl99 9). but the bulk 
of the SNe la sample have peak brightnesses within a 0.5 mag. 
range in B and V. 

An interesting feature of SNe la is that the luminosity is 
set not by the explosion, for which the deposited energy 
goes into expansion, but rather by the decay of 56 Ni and 
56 Co that are formed du ring the nucleosynthesis JArnetll982t 
iPinto & Eastma nl200ClL At the time of peak luminosity, 56 Co 
has not yet decayed, and hence the peak luminosity is a mea- 
surement of the amount of 56 Ni synthesized during the ex- 
plosion. This amount presumably depends on the progeni- 
tor and on the details of the explosion. In particular, it has 
long been known that the amount of 56 Ni synthesized depends 
in part on the asymmetry between neutrons and protons in 
the progenitor white dwarf jTruran. Arnett. & Cameronll967t 
lArnett. Truran. & Wooslevl 19711) . 

In this Letter we explore how the intrinsic variation in the 
initial CNO abundances of the progenitors translates into a 
variation in the mass of 56 Ni synthesized (M( 56 Ni)) and hence 
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in the peak luminosity of SNe la. We show analytically, using 
mass and charge conservation, that M( 56 Ni) depends linearly 
on the initial metallicity Z, and that the observed variation 
in the Z of field dwarfs, from I to 3 times solar, leads to a 
« 25% variation in M( 56 Ni), with most of the effect occurring 
at Z > Zq . This is the expected variation if SNe la progenitors 
uniformly represent the range of CNO abundances observed 
in stars today. Our conclusion rests on (1) considerations of 
nuclear equilibrium during the explosion, and (2) the observed 
scatter in Z of stars within the Galactic disk. 

The dependence of the ejected 56 Ni mass on the pro- 
genitor's initial Z has been investiga ted previously, both 
via the evol uti on of the progenitor ilUmeda et alJ 1 1999a: 
l Umeda et alJ Il999bt iDomfnguez. Hoflich. & Straniero 
1200 11 hereafter DPI) and the ex plosion itself 
(Hofli ch. Wheeler. & Thielemann! 119981 llwamoto etafl 
1999; D01). Although these one-dimensional simulations 
are sophisticated in their treatment of the flame microphysics 
and the nuclear burning, we desire to elucidate features 
that are robust to any complicated hydrodynamics. Our 
demonstration that the mass of nickel depends linearly on the 
initial Z of the progenitor serves as a check and stimulus for 
future numerical studies and provides insight into possible 
evolutionary effects. 

In §|3 we demonstrate the linear dependence of M( 56 Ni) on 
the progenitor's Z, and compare this with detailed nucleosyn- 
thesis calculations and previous numerical work. We spec- 
ulate on the source of differences between different calcu- 
lations, and then discuss the intrinsic scatter in CNO abun- 
dances in the ISM (§01 and the implications (§f4jl of our re- 
sults. 



2. VARIATIONS IN THE MASS OF NICKEL EJECTED 
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Nearly all one-dimensional Chandrasekhar mass models of 
SNe la produce most of their 56 Ni in a burn to nuclear statisti- 
cal equi librium (NSE) between the mass shel l s 0.2 and 
0.8M ? JNomoto. Thielemann. & Yokoill984t iHtiflich et al.l 
1 199a llwamoto et alJll999t IHtiflich et al.ll2000l). In this re - 
gion, unhke in the innermost 0.2M Q ( Brach witz et al]l2000ft . 
weak interactions operate on timescales longer than the time 
for the thermonuclear burning front to disrupt the white dwarf. 
Following this rapid burn to NSE, most of the mass is in the 
iron-peak nuclei 56 Ni, 58 Ni, and 54 Fe. First consider the case 
when 56 Ni and 58 Ni are the only two competing species. The 
addition of 54 Fe will be considered next. Mass and charge 



conservation, 



i=l 



(1) 



imply that the mass fraction of 56 Ni depends linearly on Y e , 

X( 56 Ni) = 587,-28, (2) 

where isotope i has Z,- protons, A,- nucleons (protons + neu- 
trons), and a mass fraction Xj. The aggregate ensemble has a 
proton to nucleon ratio of Y e . 

Most of a main-sequence star's initial Z comes from the 
CNO and 56 Fe nuclei inherited from its ambient interstellar 
medium. The slowest step in the hydrogen burning CNO cy- 
cle is proton capture onto 14 N. This results in all the CNO cat- 
alysts piling up into 14 N when hydrogen burning on the main 
sequence is completed. During helium burning the reactions 
14 N(a,7) 18 F(/3V e ) ls O(a,7) 22 Ne convert all of the 14 N into 
22 Ne. Thus, the mass fraction of 22 Ne in the carbon-oxygen 
white dwarf remnant is 
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where the mass fractions refer to the original distribution of 
the star (prior to main-sequence burning). For a uniform dis- 
tribution of 22 Ne and 56 Fe throughout the star, eq. ([0 gives 
the initial Y e of the white dwarf 
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According to the argument at the beginning of this section, 
this Y e is fixed in the region where the 56 Ni is created. Substi- 
tuting eqs. @ and (0 into eq. (0 leads to a linear expression 
for the mass fraction of 56 Ni in an NSE distribution in terms 
of the main-sequence star's initial Z 
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X( 56 Ni)=l-58 



12 



14 



16 



'Q) + XC 



28 



= 1 



Z 

-0.057—. 

z 



(5) 



The average peak B and V m agnitudes of nearby SNe la 
(Gibs on et aflboOOt ISaha et alJll999h strongly imply that a 
fiducial SNe la produces sa 0.6M Q of 56 Ni. Taking eq. l|5} 
to represent the mass fraction of 56 Ni relative to this fiducial 
mass gives 



M( 56 Ni) w O.6M ( 1-0.057 



(6) 



Here we assume that Y e is uniform throughout the star, and 
that all material within the 56 Ni-producing mass shell passes 
through NSE with normal freeze-out. 



If 54 Fe is also present then an add itional Saha-l ike equa- 
tion for the chemical potentials ( Cliff ord" & Tavleri n"965) is 
required. Consider the chemical equation a 56 Ni + /3 54 Fe <-> 
7 58 Ni. Both 54 Fe and 58 Ni carry two extra neutrons; this 
fixes the stoichiometric ratios 7//? = 1 and (3/a= 14. As a 
result, the equation for chemical equilibrium is of the form 
X( 54 Fe)/X( 58 Ni) = [/(r,p)/X( 56 Ni)] 1 / 14 . Here f(T,p) con- 
tains the mass excesses and phase space factors common to 
the chemical potential of the three species. Expanding about 
T = 5 x 10 9 K, p = 10 8 gem" 3 , and X( 56 Ni) = 1 .0, we find that 
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For this ratio, X( 54 Fe)/X( 58 Ni) » 2, the relationship be- 
tween Z and X( 56 Ni) becomes slightly shallower than eq. 
X( 56 Ni) = l -0.054(Z/Z Q ). Similai- considerations show that 
56 Fe is not present unless Y e is much smaller than the case 
considered. Equation Q does not hold if an alpha-rich freeze- 
out occurs, in which case (a, 7) reactions convert 54 Fe to 
58 Ni. Hence the ratio X( 58 Ni)/X( 54 Fe) is sensitive to the flame 
speed. The 56 Ni mass is affected only slightly, however, as 
the charge-to-mass of 58 Ni (which sets the slope in eq. Q) is 
28/58 = 0.483 and differs by only 0.3% from the charge-to- 
mass of 54 Fe, 26 /54 = 0.48 1 . Our main result, eq. (0, is thus 
relatively insensitive to the precise details of the explosion, so 
long as NSE favors a 56 Ni-dominated "iron" peak. 

The simple linear relation between Y e and the mass fraction 
of 56 Ni in eq. Q is robust, as it only relies on basic properties 
of NSE. As a test of this analytical resu lt, we cal culated the 
mass of 56 Ni ejected by W7-like models jNomoto et alll984t 
iThielemann. Nomoto. &Yokoil986|: llwamoto et all 1 9 99) by 
integrating a 510 isotope nuclear reaction network over the 
thermodynamic al trajectories for a wide range of Z. The re- 
sults are shown in Fig. ^ with the short-dashed curve. We 
also plot the linear relation (eq. |6); solid curve), with Z/Z 
adjusted to give the 22 Ne abundance used by W7, X( 22 Ne) = 
0.025Z/Zq. Most of the differences between eq. and the 
detailed W7-like models is attributable to our assumption that 
all of the 56 Ni comes from the 0.2 to 0.8M© region in the 
white dwarf, with an additional small correction for our ne- 
glect of weak interactions that slightly decrease Y e . The differ- 
ence between the slope given by the detailed W7-like models 
and the analytical model is less than 5%. 

Previous investigations ( Htiflich et a DEM! D0l) into the 
effect of a varying X( 22 Ne) do not agree on the mass of 56 Ni 
produced. Htiflich et al. ( 1998) found that varying the metal- 
licity from Z = 0.1Z Q to Z = 10Z Q produced only an « 4% 
variation in the 56 Ni mass ejected, in sharp contr ast to our re- 
sults as well the results of Iwamoto et al. ( 1999), who found 
that increasing Z from zero to solar in fast-deflagration mod- 
els leads to an w 10% decrease in the 56 Ni mass ejected. Both 
of these works relied on post-processing the thermodynamical 
trajectories, as we have done. Recently, D01 investigated the 
range 1O- 1O Z to 0.02Z Q . Of the one-dimensional calcula- 
tions just mentioned, only this one accounts for the effect of Y e 
on the energy generation rate. The long-dashed curve in Fig.^ 
shows the D01 results for their 1.5M progenitors; other pro- 
genitor masses in their survey display the same trend with 
Z. Our analytical result and post-processed W7-like models 
essentially agree with the findings of D01 over the range of 
metallicities common to all three calculations. As is evident 
in Fig. ^ the largest variation in the mass of 56 Ni occurs at 
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FIG. 1. — Mass of 56 Ni ejected by SNe la as a function of the initial 
metallicity Z. Shown is the linear relation (solid curve; the curvature is 
from the logarithmic abscissa) of eq. |6| for X( 22 Ne) = 0.024(Z/Zq), a se- 
quence of W7-like models (short-dashed cun>e), and the calculation of D01 
for 1. 5 Mq progenitors (long-dashed curve). Other progenitor masses in the 
D01 survey display the same trend with Z. As indicated by the arrows, a scat- 
ter of 3 about the mean in Z of the main-sequence stars that produce white 
dwarfs leads to a variation of about 25% (0. 13 Mq) of 56 Ni ejected if the met- 
als are uniformly distributed within the white dwarf. A factor of seven scatter 
about the mean in the initial metallicity corresponds to a factor of 2 variation 
in M( 56 Ni). 



Z > Zq. We note that the calculations of D01 evolve a main- 
sequence star into the SN e la progenitor, whereas our calcu- 
lation and those of Iwamoto et al. start from a given 
white dwarf configuration. 

As a caveat, we note that our post-processing of the W7 
thermodynamic trajectories is not completely self-consistent. 
The reason is that the temperature and density profiles of 
the W7 were calculated using the energy released by burn- 
ing matter of solar Z. While Z is likely to influence the 
flame propagati on (via the change in the rate of energy pro- 
duction; Hix & Thielemann 1996, 1999), in the mass range 
under consideration Si-burnin g is complete, so our assump- 
tion of NSE still holds. Indeed, Dominguez & Hoflichl ( 2000) 
found that the mass of 56 Ni synthesized was also rather insen- 
sitive to the details of the flame microphysics. It is also pos- 
sible the density, p (r , at which a transition from deflagration 
to deto nation occurs will influence the amou nt of 56 Ni pro- 
duced JHoflich. Khokhlov. & Wheeled ll 995. and references 
therein). Fits to observations seem to require, however, that 
p tr w 10 7 gcm" 3 which corresponds to a Lagrangian mass co- 
ordinate of ~ 1.0M Q , which is exterior to the main 56 Ni- 
producing layers. Further numerical studies should elucidate 
the source of the discrepancy. In particular they should de- 
termine if other effects, such as a systematic change in the 
mass coordinates where most of the 56 Ni is produced, offsets 
or enhances the variation in M( 56 Ni) that we find. The demon- 
strated linear dependence between M( 56 Ni) and Z will make 
it easier to untangle competing physical effects in a fully self- 
consistent calculation. 

3. SCATTER IN THE INITIAL METALLICITY AND THE INDUCED 
BRIGHTNESS VARIATIONS 

Stellar abundance determinations are discussed in terms 
of an elemental abundance relative to iron, [X/Fe], as a 
function of the iron to hydrogen ratio [Fe/H], primarily be- 
cause [Fe/H] is relatively easy to measure in stars. The 
[Fe/H] ratio represents a chronometer in that the accumu- 



lation of iron in the I SM increases monotonically with time 
JWheeler. Sneden. & Trurai3IT989t) . Calibration of [Fe/H] as 
a function of time forms the basis of the age-metallicity rela- 
tionship. 

The Z of local field stars rapidly increased about 10-13 Gyr 
ago during formation of the Galaxy's disk and then increased 
much more gradually over the last ~ 10 Gyr jTwarog |[1980|; 
lEdvardsson et alJl993tlFeltzing. Ho lmberg. & HurlevfeOOl . 
More importantly for our purposes, however, is the rela- 
tively large scatte r in stellar metallicit ies, A[Fe/H] ~ 0.5, 
at any given age. iFeltzing et alJ J2001I) constructed an age- 
metallicity diagram for 5828 dwarf and sub-dwarf stars from 
the Hipparcos Catalog using evolutionary tracks to derive 
ages and Stromgren photometry to derive metallicities. They 
concluded that the age-metallicity diagram is well-populated 
at all ages, that old, metal-rich stars do exist, and that the scat- 
ter in Z at any given age is larger than the observational errors - 
Other surveys of stellar metallicities (Edvar dsson et alJll993t 
IChen et all2000l) are in good agreement with these trends. 

The most abundant elements in the Galaxy after H 
and He are CNO. Both [C/Fe] and [N/Fe] in halo and 
disk dwarfs are observed to be r oughly solar and con- 
stant GEH LEU iCarbon et alJ fiw%: iWheeler et alJ 
( 1989)). The [O/Fe] ratio is larger at low metallicities — 
oxygen being the dominant element ejected by SNe II — 
and the n slowly dec re ases due to variations in mass 
and Z llGrat.tonl 1 19851 iPeterson Kurucz. & Carnevl 1199(1 
Timmes, Woosl ev! & Weaverl Il995t) . Within these general 
trends is a relatively large scatter, A [C/Fe] ~ A [N/Fe] ~ 
A [O/Fe] ~ 0.5 dex, at any given [Fe/H]. 

According to the simple analytical relation (eq. Q) and the 
detailed W7-like models (Fig. 0, a scatter of a factor of 3 
about the mean in the initial metallicity (|Z < Z < 3Z Q ) 
leads to a variation of about 25% (0.13Mq) in the mass 
of 56 Ni ejected by SNe la if the 22 Ne and 56 Fe are uni- 
formly distributed within the white dwarf. The minimum 
peak brightness variatio ns caused by this varia tion in 56 Ni 
mass are \AM V \ w 0.2 JPinto & Eastman! 120011) . Thus, the 
amplitude of this effect cannot account for all of the ob- 
served variation in peak luminosity of local SNe la (0.5 mag- 
nitude in B and V). The observed scatter in peak brightnesses 
may be even larger, as Cepheid-based distances to the host 
galaxies of peculiar events such as sub-luminous SN 1991bg 
or brighter-than -normal SN 1991T haven't been measured 
vet JSahaetalJI 19991 IContardo. Le ibundgut . & Vaccall2000t 
Leibundgut 2000). There is evidence for a larger scatter when 
more distant su pernovae are included (Hamuv et alJ 1199a 
iRiess et aflll998l) . It would take a scatter of about a factor 
of seven (A[Fe/H] ~ 0.8dex) in the initial Z to account for a 
factor of two variation in the 56 Ni mass and peak luminosity. 

4. IMPLICATIONS AND SUMMARY 

Using the properties of NSE, we find that where weak in- 
teractions are unimportant the mass fraction of 56 Ni produced 
depends linearly upon the initial metallicity of the white dwarf 
progenitor. This result is robust: so long as the region that 
reaches NSE does so on a timescale over which Y e is nearly 
constant, then the mass of 56 Ni produced is largely indepen- 
dent of the detailed physics of the flame front propagation. 
This fact offers a check on sophisticated numerical calcula- 
tions and could be exploited to understand better the disagree- 
ments between different codes. 

The variation from i to 3 times solar metallicity observed 
in dwarf stars in the Galactic disk implies that SNe la should 
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have a minimum variation of 25% in the ejected Ni mass. 
This has implications for the brightness variations observed 
in both near and distant SNe la and for galactic chemical evo- 
lution. The most distant SNe la observed today have redshifts 
z < 1 . This corresponds to a lookback time of about 4-7 Gyr 
depending on the cosmological model, and implies a mean 
[Fe/H] ratio between -0.1 and -0.3. There is, however, still 
a scatter of A[Fe/H] ~ 0.5 at these mean [Fe/H] ratios (see 
the discussion in § 3). We therefore expect that the variation 
in peak luminosities of SNe la in spiral galaxies at z < 1 will 
show the same minimum variation of |AAfy| ~ 0.2 in peak 
luminosity as nearby SNe la. This variation is superposed on 
other evolution ary effects, suc h as that from a reduction in the 
C/O ratio llHoflich et all 19981 D01). 

For combustion to NSE with no freeze-out (see § 2), the 
56 Ni mass is anticorrelated with the mass of 54 Fe and 58 Ni 
ejected, and subluminous SNe la will tend, therefore, to have 
larger 54 Fe/ 56 Fe ratios than brighter ones. Because SNe la 
produce > i of the galactic iron-peak nuclei 3 , the isotopic ra- 
tios among the iron group in SNe la ejecta shoul d not exceed 
the solar ratios by ab out a factor of two dlwamoto et alii 19991 
IWheeler et aflll989l) . There are several uncertainties with the 
54 Fe/ 56 Fe ratio. First, some of the 54 Fe is produced in the 
core where weak interactions are important for determining 
the final Y e (although this may be alleviated by the overturning 
of matter in the cor e from Rayleigh-Taylor instabilities; see 
iGamezo et al.l2003l for a recent calculation). At high densities 
these reaction rates are sensitive to the input nuclear p hysics 
(Brachwitz et al. 2000; Martfnez-Pinedo. Langanke. & Dean 
2000). Isotop ic measurements of co s mic rays, suc h as 
with Ulysses dConnell & Simpsonlll997l IConnellll200ll) and 
ACE iWiedenbeck et al.l 120011) can probe the evolution of 
54 Fe/ 56 Fe over the past « 5 Gyr. 



Since most of the 56 Ni created in one-dimensional SNe la 
models lies in the 0.2M Q -0.8M Q mass shell, our as- 
sumption of a uniform 22 Ne distribution deserves close 
scrutiny. There has long been speculation about sedi- 
mentation of 22 N e and its effect on the cooling of iso- 
lated white dwar fs (Bilds ten & Hall200ltlDelove & Bi ldsten 
120021 [Hansen et alJl2002l and references thereinl. If, for ex- 
ample, all of the 22 Ne from the outermost 0.6 M Q were to set- 
tle into the shell between 0.2M Q and 0.8 M b, then the effec- 
tive Y„ in this shell would double. Indeed, Bi ldsten & Halll 
d2001l> noted that the production of 54 Fe was an indirect test 
of the sedimentation of 22 Ne. For nearby SNe la, which pre- 
sumably sample a range of progenitor ages, the variability in 
56 Ni production would increase to 50%, or about AM V « 0.3. 
Since it takes about 7 Gyr for 22 Ne at the surface to fall 
through the outer sa .4M^ of a 1.2M Q CO white dwarf 
(Bildsten & Hall 2001), the effect of sedimentation will be di- 
minished for those SNe la at z ~ 1 . 
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3 SNe la produce « 0.8 Mq of iron-group nuclei while Type II SNe 
produce »s O.IMq, and the ratio of thermonuclea r to core collapse 
events is about 0.15-0.27 in the Galaxy ivan den Berah & Tammann 1991; 
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